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Ballistic transport of helical edge modes in two-dimensional topological insulators is protected 
by time-reversal symmetry. Recently it was pointed out [T] that coupling of non-interacting helical 
electrons to an array of randomly anisotropic Kondo impurities can lead to a spontaneous breaking 
of the symmetry and, thus, can remove this protection. We have analyzed effects of the inter¬ 
action between the electrons using a combination of the functional and the Abelian bosonization 
approaches. The suppression of the ballistic transport turns out to be robust in a broad range of 
the interaction strength. We have evaluated the renormalization of the localization length and have 
found that, for strong interaction, it is substantial. We have identified various regimes of the dc 
transport and discussed its temperature and sample size dependencies in each of the regimes. 

PACS numbers: 71.10.Pm, 73.43.-f, 72.15.Nj, 73.20.F 


Introduction: Electron transport in time-reversal in¬ 
variant topological insulators (TI) has become a hot topic 
of research during past several years, see Refs. [ng for re¬ 
views. The bulk electron states in the TIs are gapped, 
nevertheless, dc transport is possible since it is provided 
by low-dimensional helical edge modes. Helicity means 
the lock-in relation between electron spin and momen¬ 
tum: helical electrons propagating in opposite directions 
have opposite spins [SIT]. An elastic backscattering of 
a helical electron must be accompanied by a spin-flip. 
Therefore, the helical electrons are immune to effects of 
potential disorder such as localization. 

Recent experimental studies of the charge transport 
through ID channels at the helical edges of 2D TIs 0 - 
[TT| made of quantum-wells [121 dg demonstrated that 
the transport is indeed close to be ballistic as long as 
the samples are small. li El. However, longer edges 
exhibit lower conductance iHiinKig which is evidence for 
back-scattering. Moreover, the absence of clear temper¬ 
ature dependence of the sub-ballistic conductance in a 
broad temperature interval m suggests that this back- 
scattering is probably elastic. 

Robustness of the ballistic transport in the Tis was 
discussed in several theoretical papers. Spin-flips needed 
for the backscattering could be, for example, due to spin- 
1/2 Kondo impurities. However, the Kondo-screening 
of the spin would recover the ballistic transport at low 
temperatures dZllIS]. The inelastic processes caused by 
interactions in the presence of disorder are predicted to 
result in temperature-dependent contributions to the dc 
conductance and conductance fluctuations [TMSg . Such 
corrections are also frozen out at cooling due to the lack 
of the phase space. Thus, the explanation of the dc trans¬ 
port in the Tis, which (i) is non-ballistic at low tempera¬ 


tures, and (ii) can be temperature independent, remains 
a theoretical challenge. 

The first step toward the understanding of the non- 
ballistic temperature-independent transport through ID 
helical edges was taken in Ref. [1] where the idea of spon¬ 
taneous breaking of time-reversal symmetry was pro¬ 
posed. It was demonstrated that helical ID electrons, 
which do not interact with each other, are localized if 
they are coupled to an array of the Kondo impurities. 
This coupling can be described by the Hamiltonian: 


Hb = dxpsJ± 


{S+ + h.c. . (1) 


Here (■0- ) describes spin up right moving (spin down 
left moving) in a;-direction helical fermions 'ipR,^ 
kp is the Fermi momentum; ps is the impurity density; 
J_L = {Jx + Jy)/‘2‘ is the coupling constants between the 
Kondo- and the fermion- spins; e = {Jx — Jy)/J±. is the 
dimensionless parameter of the anisotropy in the plane of 
the2DTI (XY-plane); S'± are the Kondo spin operators. 

For isotropic couplings, e = 0, the indirect interaction 
between spins induces a slowly varying in space and time 
spin polarization. Homogeneous time-independent polar¬ 
ization would create a gap, Ag = spsJ±, in the spectrum 
of fermions (s = 1/2 is the impurity spin). Fluctuations 
of the polarization result in heavy but gapless “polaronic” 
complexes of helical electrons dressed by slow spinons. 
These complexes are charged and can support ballistic 
transport with a strongly reduced Drude weight. A ran¬ 
dom anisotropy, e{x) ^ 0, quenches the local spin polar¬ 
ization and thus causes spontaneous symmetry breaking. 
The polaronic complexes loose their protection from the 
back-scattering and undergo the Anderson localization 
with maybe large but finite localization radius, 
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FIG. 1. Typical dependence of the effective gap on the 
Luttinger parameter, see Eq. ||^. 


In this Letter, we explore the charge localization in the 
system of interacting helical electrons - Helical Luttinger 
Liquid (HLL) - coupled to the array of the Kondo im¬ 
purities, and analyze the temperature dependence of the 
dc transport. Interactions in the HLL are characterized 
by the Luttinger parameter K = {1 + (if is the 

dimensionless interaction strength). We prove that the 
moderate attraction, 1 < K < 2, and (almost) arbitrary 
repulsion, K < 1, do not change the qualitative picture 
of the non-interacting system: the effective theory, which 
describes localization in the HLL coupled to the Kondo 
array, remains valid, though the gap. A, and the local¬ 
ization radius, Lioc, are substantially renormalized: 


A(K) / Eb Lioc 

Ao ^ [k^AoJ ’ r(o) 


Aq \ 
KA{K)) 


4 

3 


( 2 ) 


see Fig{l] Here Eb is the UV energy cut-off which is of 
the order of the bulk gap in the Tl. 


The energy scales which govern different temperature 
regimes of the dc transport are sketched on Figj^ They 
are the temperature of many-body localization transition 
|24j . Tmbl, and the de-pinning energy, Epin, defined be¬ 
low, see Eq. a- We will assume that [see the discussion 
after Eq.(211] 


Tmbl Epin ^ A . (3) 


If T < Tmbl all excitations are localized and ballis¬ 
tic transport is suppressed. If Tmbl < T < Ep-^, the 
dc conductivity is finite albeit low and is of a quantum 
nature. The transport becomes thermally activated as 
T —>■ Tpin and turns into semiclassical one in the inter¬ 
val Tpin < T < A. Power-low dependencies of a^c at 
T 3> Tmbl result from the interaction-dependent renor¬ 
malization of J_L, see Eqs.(22|^|. 


The model and calculations: The Hamiltonian of the 
HLL coupled to the array of the Kondo impurities is: 
H = Ho + Hint + Hb, where the first two terms describe 
the free fermions and the interaction between them, re¬ 


spectively: 



Hint = ^ J dx {p++ p-f, P± = V'±'0±; (5) 


here vp is the Fermi velocity and v is the density of 
states in the HLL. The backscattering term Hi, (caused 
by the coupling of the electrons to the Kondo impu¬ 
rities) is defined in Eq. 0. We neglect the forward¬ 
scattering term ~ Jz^z since a unitary transformation 
of the Hamiltonian allows one to map the model with 
the parameters {K, ^ 0} to its counterpart with the 

effective Luttinger parameter K = K{1 — J^vI2K')'^ and 
Jz = d [25l [26] . Thus, Hint takes into account both the 
direct electron-electron interaction and the interaction 
mediated by the z-coupling to the Kondo impurities. 

Let us briefly recapitulate important steps, equations 
and basic assumptions from the paper [Tj. The spin or¬ 
dering allows us to develop a path integral formulation 
of the problem using the parametrization of each spin by 
its azimuthal angle, a, and projection on z-axis, n^: 

S^{x,t) = - nl{x,T ), jn^j < 1; (6) 

T denotes the imaginary time. It is convenient to use 
the spins rotation: 5'^(a:,—)■ S^{x,t) which 
leads to the redefinition of the anisotropy parameter: 
^^ 2 ikpx pgj. high density of the Kondo 

array with a weak irregularity in positions of the spin 
impurities, the model can be simplihed by treating 
and as random variables: real Ps(x) and com¬ 

plex e{x). Fluctuations of ps do not effect the dc trans¬ 
port in the HLL, cf. Ref.[T], thus, one can use the aver¬ 
aged density. On the contrary, randomness of e{x) plays 
the crucial role. Without loss of generality, e(a;) can be 
treated as a Gaussian random function with zero mean 



FIG. 2. (color on-line) Sketch of the temperature dependence 
of the dc conductivity in different regimes, see explanations 
in Sect. “TG transport” and Eqs. ( |22p3j l in the text. 
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value and short-range correlations 


(e(a:))dis = 0; (e(a:)e*(a;'))dis = —6{x - x'). (7) 

Ps 


As a result, the Lagrangian density for the HLL coupled 
to the Kondo array can be presented as: 


C = 'tp^ 


d+ J 
J* d- 


J = 


nt 


(e— + ee*“) ; 


ispsU^dra-, ( 8 ) 


where d± = =p ivpd^ are the chiral derivatives; 

= {V’+, is the fermionic spinor held. Further cal¬ 
culations are based on the scale separation: the fermionic 
variables are much faster than the spin ones. Therefore, 
as will be verihed below. 


< 1, A ^ = (9) 


To describe the composite fermion-spinon excitations, 
we perform a gauge transformation of the fermionic helds 
—>■ ipri- After a point-split of the interaction 
term and evaluation of the Jacobian, the Lagrangian den¬ 
sity acquires the form: 

C[p),a]~ C\j^jeic A -; ( 10 ) 

where v = vp/K is the renormalized velocity, and sub¬ 
leading terms Prjdr^a are neglected in L due to inequal¬ 
ities Eq.(|^. 

We start our analysis of the effects of the electron- 
electron interaction with the zero anisotropy case, e = 0. 
The classical conhguration of the spin variables is 

= 0, = const. 


This expression has been obtained after neglecting small 
and slow spatial huctuations in Ps{x) and nz{x). Be¬ 
sides, Eq. 0 is valid only provided that 

T < A{S) <£<^Eb. 


The renormalization stops as soon as the energy £ be¬ 
comes smaller than A since the fermions are massive. 
Accordingly, the effective fermionic gap Af can be de¬ 
termined from the self-consistency equation: 


Af = A(£) 


e=Ai 


( 12 ) 


Equation (12) leads to the result Eq.([^ for attractive, 
1 < iF < 2, or weakly repulsive, 0 < 1 — AT ^ 1 in¬ 
teractions. However, this straightforward way of calcula¬ 
tions suffers from a serious drawback since the renormal¬ 
ization procedure dynamically generates potentially rel¬ 
evant multiparticle scatterings. It is technically difficult 
to analyze their relevance within the functional bosoniza- 
tion approach. To justify the self-consistent derivation of 
Af and to analyze the case of stronger repulsion, one 
needs an alternative approach. 

Such an approach can be based on the bosonization of 
the fermionic part of C followed by the analysis of the 
fully bosonized Lagrangian: 


C[p, a] ~ £sg + 


V jdxa) 
8 ttK 


- ispsTizdra ; (13) 


where at e = 0 


£-SG = K,v)-\ - ^ -cos(2^); (14) 


2'Ka 




1 

2'kKv 




■ (15) 


Our goal is to derive an effective theory describing fluc¬ 
tuations around the classical solution. For g = 0, this 
is achieved by integrating out all massive modes: firstly, 
the fermions with the gap Aq, and, secondly, the non- 
dynamical variable (in the quadratic approximation). 
The interaction can be taken into account with the help of 
the functional bosonization, which involves the Hubbard- 
Stratonovich decoupling of the interaction term and the 
gauge transformation of the fermionic fields [23 [28] . Fol¬ 
lowing these standard steps, one can show that the main 
non-perturbative effect of the weak interaction is the 
renormalization of the backscattering amplitude, which 
acquires an effective energy-dependence, J±{£). Similar 
to the bare relation Aq (x Jl , we introduce the renormal¬ 
ized gap at each step of the renormalization procedure: 
A oc J±{£). Thus, the gap also becomes energy depen¬ 
dent. The energy-dependence reads as: 


Here a ^ vp/Ep is the smallest spatial scale (the UV 
cut-off); (/) is a composite phase: p = p — a/2, with p 
being the usual bosonic phase - its gradient is related to 
the fermionic density: dx4> =—'^{p++ P-) E3- 


Similar to Eq.(lO), sub-leading gradients of a are ne¬ 


glected in Eq.(13). This can be justified since £sG cor¬ 


responds to the quantum Sine-Gordon model where the 
relevant vertex ~ Aq cos(2(/) generates the bosonic mass 
at AT < 2 [30] resulting in the scale separation between 
bosonic and spin degrees of freedom. The effective gap 
A in the spectrum of bosons can be determined with 
the help of the Feynmann variational method [22] • After 
neglecting fluctuations of Uz, the gap equation can be 
written as: 


Eb J Eb V J 


K 


(16) 


A(£:) = Ao X 



( 11 ) 


Equation (16) is valid at T < A for AT < 2, i.e., 
for arbitrary strong repulsion as well as for weak and 






















4 


moderate attraction m- The multiparticle backscat- 
tering, which is described by (less relevant) higher ver¬ 


tices ~ cos 


2n(j) 


in £sg, n- = 2,3,..., can be 


included into the gap equation. Such vertices yield 
corrections to the RHS of Eq.(16l of order 0{'n? x 


/Eb)x{KA/Eb)'^ ^). Assuming that Eb 

and A Eb, contributions from the multiparticle 
backscattering are sub-leading and can be neglected at 
K>l/logiEB/Ao) [32]. 

At AT ~ 1, Eqs.( f2) and (16) coincide - the bosonic gap 


A and the gap of the interacting helical fermions Af are 
equivalent for weak and moderate interaction strengths. 
The fluctuations of Uz result in small fluctuations of A. 
One can analyze a change of the free energy caused by the 
gap fluctuations and prove an equivalence of the bosonic- 
and fermionic gaps for arbitrary K [33] . 

The effective low energy theory for the phase a can 


now be derived either from the Lagrangian (10) or from 
its fully bosonized counterpart We use the former 

option, namely, we restore the anisotropy, neglect g in 


the hrst term C of (10) but take into account the inter¬ 


action induced renormalizations, i.e., we use the effective 
gap A, the effective velocity v, and the renormalized 
factor vISttK in the gradient term rather than the bare 
values. We will use A in further calculations since, un¬ 
like Af, it allows us to describe the strong repulsion. 

Integrating out the massive variables and Uz, 

we obtain: 


Eql -^ll , "Uq, ) 

V -1- c.c.) ; 

(17) 

where we have defined 

I A^ / 

Eb\ 

(18) 

V = -bg 

47ra Eb \ 

—j, aV < A; 

aD 

V AK Ay/K 

< I. 

(19) 


At K = 1, Eqs.(17jl9|) are equivalent to the results of 


Ref.[T]. The composite particles are slow in the absence 
of the interactions, Va{K = l)/nF ^ 1- One can see 
that, as K decreases, the ratio Va/v increases remaining 
small as long as AT > 1/log(AF/Ao). Inequality v ^ 
Va reflects the scale separation between the interacting 
(fast) helical fermions and the (slow) composite particles. 
Eq.(17l is valid at AT < 2, therefore ATq, ^ aD^fKjA is 
small with- and without interactions. 

DC transport: A regular anisotropy, e(x) = const, 
would pin the phase, a ~ arg(e). The pining is similar 
to the effect of a magnetic field applied in the XY-plane: 
it breaks time reversal symmetry and opens a global gap 
in the spectrum of the composite quasiparticles trivially 
blocking the ballistic dc transport. 

Random fluctuations of arg(e), Eq. 0. prevent the 
global gap from opening but are able to localize the com¬ 
posite particles at T —>• 0. Indeed, Eq. 0 describes 


a disordered Sine-Gordon model for which the localiza¬ 
tion was demonstrated long ago [34]. Since <C 1, 
the localization length can be evaluated by the standard 
optimization procedure [55113S] : 


(») 

w is defined in Eq.Q. Let us now discuss different tem¬ 
perature regimes of the dc transport. 

The localization strongly effects the low temperature 
dc transport in not too short samples, L > Lioc: the con¬ 
ductance is resistive (finite but smaller than the ballistic 
one) and temperature-independent in transient samples 
with h ^ Aloe, and the conductance vanishes in long 
samples, A ^ Aioc. 

Sizable dc transport in long samples (A ^ Aioc) ap¬ 
pears only at higher temperatures: if Tmbl <T < Apin, 
the fermions are still gapped but many-body states be¬ 
come delocalized and are able to support weak dc quan¬ 
tum transport similar to transport in glassy systems. As 
T —>• Apin, the transport becomes classical. The straight¬ 
forward estimate yields (cf. Ref. [55] 1: 

)“''“. ( 21 ) 


The classical de-pinning energy is supposed to exceed the 
quantum energy scale, Ambl ^ Apin (a rigorous theory 
for Tmbl at ATq, ^ 1 is missing). Semiclassical con¬ 
ductivity of the composite particles can be estimated as 
CTc OC [36], where ~ (T/Epin)'^{KaEpin)~^ 

is an effective temperature dependent transport time m- 
Using Eq.(19), we arrive at 


CTc OC 


VF 

K K ■ 


K 


T 

Er\\x 


( 22 ) 


The factor (ATc/AT)^ «C 1 reflects the suppressed Drude 
weight. Equation ( [^ is valid provided that Apin < T «C 

A. 


Strictly speaking, Eq.(16l is valid at T < A. Nev¬ 
ertheless, we can draw some qualitative conclusions for 
T ~ A ^ Eb, where the dc conductivity is dominated 
by thermally activated fermionic quasiparticles and can 
be estimated as 1551: 


CTf OC - vpTo ; (23) 

where tq is governed by the disorder of the Kondo lat¬ 
tice, i.e., by the randomness in p^. The theory for tq is 
beyond the scope of the present Letter. If A ^ T ^ Eb, 
the gap becomes temperature dependent and shrinks. As 
a result, our theory looses its validity and the dc trans¬ 
port should reflect different physics. 
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Validity: Our consideration is based on several assump¬ 
tions: Firstly, the Kondo array is dense, the bare coupling 
constant is small and the XY-anisotropy is weak 


Psa ^ 


1 , 


i'J± <C 1, |e|, re <C 1. 


(24) 


Combining these inequalities with Eq.Q, one can check 
that A ^ Eb and justify Eq.(|^. 

Secondly, we have neglected the Kondo effect which 
is permissible only provided that A exceeds the Kondo 
temperature of a single Kondo impurity embedded into 
the HLL, Tfc- The standard estimate for obtained 
from the RG approach reveals two regimes: ~ 

Eb exp{—l/vJ±) is exponentially small at K —>■ 1 but 
becomes larger due to the renormalization of J± in the 
interacting case, ^ EBi^J ^/at 

1 — K ^ vJb [13 [3H]- Comparing Tk with A from 
Eq.([^, we find that Tk is always the smallest scale. 
Therefore, for a dense Kondo array, the Kondo screening 
can be neglected. 

Finally, taking into account restrictions on K dis¬ 


cussed after Eq.(16), we determine the validity range for 


the interaction strengths: 


(log(F;B/Ao)) ^<K<2. 


(25) 


These conditions are not too restrictive, i.e., the pre¬ 
sented theory is valid in the broad range of interaction 
strength. 

Conclusions and open questions: We have demon¬ 
strated that the localization of the ID helical electrons 
coupled to the random array of the Kondo impurities 
is a generic and robust phenomenon which takes place 
in the broad range of the electron interaction strengths, 
Eq. ( |2^ . This confirms a qualitative conjecture of Ref. [T] . 
We have found and quantitatively described strong (non- 
perturbative) renormalizations of physical parameters, 
Eq.([2]). 

We have discussed possible manifestations of localiza¬ 
tion in the dc transport at low temperatures T <C A (A 
is the gap in the electron spectrum caused by local spin 
ordering), when the current is carried by slow compos¬ 
ite spinon-fermion excitations. A random anisotropy of 
the electron-spin coupling pins spin ordering and local¬ 
izes low-energy excitations. The localization length Lioc 


is given by Eq. (20). The localization is known to lead to 


an insulating state, which persists in a finite temperature 
interval T < Tmbl, where Tmbl is the temperature of 
a many-body localization transition [24] . Tmbl is ex¬ 
pected to be small compared to the classical de-pinning 


energy Epin ^ A, Eq.(21). This leads to the following 
predictions for short samples, L L\oc: the deviation of 
the conductance from the universal ballistic value should 
be small as L/L\oc- The low temperature transport of 
the samples with L ^ Lioc is characterized by strong 
mesoscopic fluctuations. Longer samples, L ^ Lioc, are 
almost prefect insulators as long as T < Tmbl- The 


theory for a temperature interval Tmbl < T < Epin is 
yet to be developed for both long and short samples. At 
higher temperatures, Tpin < T ^ A, the dc transport 
becomes semiclassical with a ^ Eq.(221, and reduced 

Drude weight {Ka/KY 1, see Eq.(19) (and Eq.(3) in 
Ref.d]). 

To develop the analytical theory, we had to restrict our 
choice of parameters to those given in Eq.(24). In par¬ 
ticular, we have assumed a small anisotropy. This choice 
results in low Tmbl and large Tmc- However, a strong 
XY-anisotropy looks more natural for the edge transport 
and it could yield stronger localization with larger Tmbl 
and smaller Tmc. Therefore, we believe that the con¬ 
sidered mechanism of suppressing the ballistic transport 
can be relevant for realistic samples. In particular, the 
temperature independent transport observed in Ref.|16j 
might be associated with the described above resistive 
regime of relatively short samples at T < Tmbl- 
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be published elsewhere. 
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Combining all together we hnd the estimate 
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